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Abstract 

We argue that coherent oscillations of the axion field excited by the misalignment mechanism 
and non-thermal leptogenesis by the saxion decay can naturally explain the observed abundance 
of dark matter and baryon asymmetry, thus providing a solution to the baryon-dark matter 
coincidence problem. The successful axionic co-genesis requires a supersymmetry breaking scale 
of 0(1O 6 ~ 7 ) GeV, which is consistent with the recently discovered standard- model like Higgs 
boson of mass about 126 GeV. Although the saxion generically decays into a pair of axions, 
their abundance sensitively depends on the saxion stabilization mechanism as well as couplings 
with the Higgs field. We discuss various ways to make the saxion dominantly decay into the 
right-handed neutrinos rather than into axions, and show that the abundance of axion dark 
radiation can be naturally as small as AiV e fj < 0(0.1), which is allowed by the Planck data. 
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I. INTRODUCTION 



There are many coincidences in nature. Some of them may be just a coincidence or 
may be a consequence of anthropic selection, while others may be due to a novel physical 
mechanism. We consider the last possibility, since it will provide us with a hint of physics 
beyond the standard model (SM). In this paper we focus on the baryon-dark matter 
coincidence problem, i.e., why the dark matter density is about five times larger than the 
baryon density, Q c ~ 5Qb- If these two have a totally different origin, it is a puzzle why 
they happen to have the density in a similar size. 

There have been many works on this long-standing puzzle. Among them, the asym- 
metric dark matter scenario has recently attracted great attention jlj . If the dark matter 
particle has asymmetry which is comparable to the baryon asymmetry, the coincidence 
problem could be solved for the dark matter mass of order the nucleon mass. Of course 
this is just a replacement of the problem. The similarity of the baryon and dark matter 
abundances is replaced with the similarity of the nucleon mass and the dark matter mass. 
So, if the dark matter mass is independent of the QCD scale, the puzzle remains unsolved. 
If the coincidence is due to a novel physical mechanism, therefore, the mechanism may be 
such that the dark matter abundance depends on the QCD scale, which determines the 
proton mass and thus the baryon density. 

The axion is a pseudo-Nambu-Goldstone boson associated with the spontaneous break- 
down of the Peccei-Quinn (PQ) symmetrvand it elegantly solves the strong CP problem 
when it couples to the QCD anomaly 2|, [3]. The invisible axion is one of the plausible 
candidates for dark matter. Interestingly, its mass and abundance depend on the QCD 
scale and the PQ breaking scale. To our knowledge, the axion is the only dark matter 
candidate whose abundance naturally depends on the QCD scale. 1 Therefore, we consider 
a possibility that the PQ sector is responsible for both the baryon asymmetry and dark 



1 Although its dependence on the QCD scale is not exactly same as the baryon density, we will see that 
it may be related to another interesting coincidence between the QCD scale and the light quark mass 
(or the weak scale). 



2 



matter, i.e., the axionic co-genesis. 

One particularly attractive mechanism that offers an explanation for the observed 
baryon asymmetry is leptogenesis ^,5]. If the right-handed neutrinos are charged under 
the PQ symmetry, the axion as well as its scalar partner, saxion, are naturally coupled 
to the right-handed neutrinos. 2 It is known that the saxion tends to dominate the energy 
density of the early universe. Then, the saxion decays into right-handed neutrinos, whose 
CP- violating decays create a lepton asymmetry 0, S]- As we shall see shortly, the saxion 
mass should be of O(l0 6 ~ 7 ) GeV, for successful non-thermal leptogenesis. 

The saxion generally decays also into a pair of axions. Such non-thermally produced 
axions remain relativistic until present, giving contributions to the effective number of 
neutrinos, N e g. 3 Some of the recent CMB observations gave a slight preference to the 



excess of N e g, coined dark radiation 
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26]. On the other hand, the Planck results did 



not confirm such excess, and place a constraint on N eS [27| 



N cS = 3.30+^ (95%; Planck + WP + highL + BAO). 



(1) 



As we shall see shortly, although the axions are generically produced by the saxion decay, 
its abundance sensitively depends on the saxion stabilization mechanism, and can be 
naturally as small as AiV e ff = 0(0.1), which is still allowed by the Planck data. 

Most of our results in this paper can be directly applied to non-supersymmetric cases, 
but we use a supersymmetric (SUSY) language for simplicity. We shall return to the 
non-SUSY case in the last section. In fact, there is an interesting implication in the 
SUSY framework. In SUSY, the saxion acquires a mass m a from the SUSY breaking. 
For a generic Kahler potential, we expect m a ~ where m 3 / 2 denotes the gravitino 
mass. Interestingly, high-scale SUSY breaking of (9(10 6-7 ) GeV is consistent with the 



2 In this case, the axion becomes a majoron 

3 The abundance of relativistic axions produced by flaton decays was studied in Ref. [9(. The late- 
time increase of N g by d ecay ing particles (e.g. saxion to two axions, and gravitino to axion and 
axino) was studied in Ref. 10(. The string axion production from the modulus decay was considered 



Refs. 
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13j | . See also Refs. l4423j and many others. The possibility that the X particle was in 



thermal equilibrium was studied in Ref. 24 1 



3 



recently discovered SM-like Higgs boson of mass about 126 GeV [28|,|29( (see Refs. [30l.l3lj). 
In fact, the observed Higgs boson mass as well as no experimental signatures of SUSY 
leads to another puzzle: if SUSY is realized in nature, why its breaking scale is much 
higher than the weak scale, which generically requires severe fine-tuning for the correct 
electroweak symmetry breaking. Our framework, the axionic co-genesis, may thus provide 
explanations for the puzzle, because the non-thermal leptogenesis by the saxion decay 
would be impossible otherwise. 

The rest of this paper is organized as follows. In Sec. Ullwe describe our set-up, and then 
discuss the axionic co-genesis of baryon, dark matter and radiation in detail in Sec. IIHI 
The saxion stabilization is discussed in Sec. IIVI The last section is devoted to discussion 
and conclusions. 

II. SET-UP 

The model we shall discuss is a PQ invariant extension of the supersymmetric SM 
where PQ charged right-handed neutrinos (Ni with i = 1,2,3) are added to implement 



the seesaw mechanism 32|. Thus N becomes massive after the PQ symmetry breaking, 
and couples to the saxion with a coupling proportional to its mass. Here and in what 
follows, we assume that the lepton number is explicitly broken by adding couplings of the 
PQ scalars with Higgs fields (such as ffl6|) or the second term of (jlOp ) in order to get rid 
of a massless majoron jf]]. 

We assume that there are multiple PQ scalar fields, X a , which acquire a vacuum 
expectation value (VEV), (X a ) = v a . The X a can be expanded around its VEV as 

X. = + (2) 

where a a and a a denote the radial and phase component of X a , respectively. The axion 
a and saxion a are expressed in terms of those components of the PQ breaking fields, 

(3) 



—f~ °a, (4) 



/ 
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where q a is the PQ charge of X a , and the PQ scale / is determined by 

/ 2 = E«- ( 5 ) 

a 

The interactions relevant to our discussion are written as 

C = -jLj a{d»a)d^a + Qa^A^ + i^AyV* + h^H^Lj + h.c.J , (6) 

below the PQ breaking scale /. Here the term generates Dirac neutrino masses after 
electroweak symmetry breaking. The saxion couples to the axions with a coupling {33) ]. 

Here one should note that both x and / depend on the charge normalization, but the 
combination x/f does not. In the following analysis, we assume a charge normalization 
such that all q a are order unity. 

The right-handed neutrinos obtain masses from the PQ breaking, which depend on 
how they couple to the PQ breaking fields. For instance, one can take the PQ charge 
assignment that allows the superpotential terms 

W = K,X a NN + k'NLH U} (8) 

where the family indices have been omitted. Then Mj ~ After integrating out the 
heavy right-handed neutrinos, light neutrino masses are generated as 

m v -, (9) 

K J 

where v is the Higgs vacuum expectation value. We are interested in the case where at 
least one of the right-handed neutrinos has a mass much smaller than the PQ breaking 
scale so that it can be produced by saxion decay. As we shall see shortly, the reference 
values are / ~ 10 12 GeV and the saxion mass m a ~ 10 6 GeV. This requires k to be 
smaller than m a /f ~ 10 -6 , which can be achieved if Aj are charged under some flavor 
symmetry. Another way to provide small masses to the right-handed neutrinos is to 
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consider higher-dimensional terms, 4 

X 2 X 

W = k—^-NN + k' —^-NLH U , (10) 

by assigning appropriate PQ charges to the involved fields. Here X a can be in general 
different from Xp, but we assume (X a ) ~ (Xp) for simplicity. In this case, the right- 
handed neutrinos obtain masses about nf 2 /A, leading to the light neutrino mass, 

k' 2 l) 2 

The value of m v is determined by v 2 /A, independent of /, in contrast to the previous 
case. The cut-off scale A is considered to be around the GUT scale (or the Planck scale 
if k *C 1), for which one can easily obtain the hierarchy M, <C / as well as small neutrino 
masses consistent with the observational bounds. Note that the axion coupling with left- 
handed leptons at low energy scales is suppressed in this case if a — (3, in contrast to the 
previous case and the original Majoron model. 

The saxion coupling to Aj is proportional to Mj, and the numerical coefficient depends 
on the model. To be concrete, we take the superpotential (|T0l) . Then the relation between 
Aj and Mi is fixed to be, 5 

A, = ^p, (12) 

where q denotes the PQ charge of the PQ breaking field that is responsible for the right- 
handed neutrino mass. From the saxion couplings x and Aj, one finds that the branching 
ratio for a — > AjAj is written 

2/2Ma 2 A /2Ma 2 \ 3/2 ^ 0.36 



Bi(a^N l N t ) = -( 1- -) B a <—B a , (13) 



4 One may instead consider the Kahler potential term, (X^/X^NN, for the PQ fields having (X a ) ~ 
{Xp). Then the right-handed neutrino masses are induced after SUSY breaking via the Giudice-Masiero 



mechanism 



34j . and naturally similar to the saxion mass. In this case, the saxion coupling for aNiNi 
is still proportional to Mi/ f with a coefficient of order unity as long as a ^ (3. 
5 In the model ([8]), is two times smaller. Our main results hold also in this case without significant 
modifications. 
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if the process is kinematically allowed. Here m a is the saxion mass, and we have defined 
r = x 2 /q 2 , which is independent of the charge normalization. B a denotes the branching 
ratio for the saxion decay into an axion pair, 

a_ r (7 647r/ 2 ' 1 } 

with T a being the total saxion decay rate. The above shows that the upper bound on 
Br(cr — > iVjiV,) is set by B a and r. If r is comparable to or smaller than 0.3, Br(cr — y 
NiNi) > B a is allowed, and the saxion can dominantly decay into right-handed neutrinos. 
In Sec. IIVt we will show that such a value of r is obtained naturally in supersymmetric 
axion models. 

Here and in what follows we do not consider the saxion decay into right-handed sneu- 
trinos, assuming that the decay is kinematically forbidden. Since the saxion mass arises 
from the SUSY breaking, it can be comparable to the right-handed sneutrino mass. Even 
if the decay mode is allowed, the following argument of non-thermal leptogenesis remains 
almost intact. 6 

Finally we mention the PQ mechanism solving the strong CP problem. In order for 
this mechanism to work, the PQ symmetry should be anomalous under the QCD so that 
the axion couples to gluons via 

c = wf°^> (15 > 

where the axion decay constant is given by F a = \[2f /A^>w- The domain wall number 
A^dw is determined by the PQ anomaly coefficient, and it counts the number of discrete 
degenerate vacua of the axion potential. For the anomalous coupling, we need either to 
assign PQ charges to the Higgs doublets or introduce additional colored matter fields 
carrying PQ charges. In the former case, an effective higgsino fi parameter is generated 



after the PQ symmetry breaking, for instance, from the superpotential 35] 



W = K H ^H u H d . (16) 



6 The LSP abundance can be affected. However, they wifl decay into the SM particies if the R-parity is 
violated. 
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It is important to note that the above superpotential induces saxion couplings to higgsinos 
as well as the SM fermions through mixing with the Higgs, which can significantly suppress 
B a depending on the values of m a and the \x parameter. This implies that B a + ^2 i Br (a — > 
NiNi) needs not to be equal to one even for the case with r of order unity. Alternatively 
one may include color-charged matter fields $ + $ c which obtain large masses from the 
coupling X Q <£>$ C in the superpotential. Unless the coupling is suppressed, the saxion decay 
into $ and $ c can be kinematically forbidden so that B a + J^B^er — > NiNi) becomes 
close to unity. 

III. COSMOLOGY 
A. Saxion dynamics 

Let us begin by discussing the cosmological role of the saxion. In the SUSY framework, 
the saxion corresponds to a flat direction in the supersymmetric limit. During inflation the 
potential receives Hubble-induced corrections, and therefore the minimum of the potential 
is generally shifted away from the true vacuum. After inflation the saxion starts to oscillate 
with a large amplitude when the Hubble parameter becomes comparable to its mass. 
Furthermore, the saxion dynamics can be significantly affected by thermal effects, and 
it is possible that the large saxion amplitude is realized by thermal corrections through 
couplings with the heavy PQ quarks [36]. Thus, the saxion tends to dominate the energy 
density of the universe unless the reheating of the inflaton is extremely low. 

In a class of the saxion stabilization models and also in the non-SUSY case, the initial 
oscillation amplitude of the saxion is considered to be of order the PQ breaking scale F a . 
The typical scale of F a in our scenario is much smaller than the Planck scale, and so, it 
is non-trivial if the saxion can dominate the universe before the decay. Even in this case, 
if the saxion sits near the origin after inflation, thermal inflation may take place and the 
saxion-dominant universe is realized. The duration of the thermal inflation does not have 
to be very long, because we do not aim at solving the moduli problem. Note however that 
in this case one needs to arrange Abw = 1 to avoid the domain wall problem, and that the 
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axion decay constant is constrained to be smaller than about 4 x 10 GeV since otherwise 



the axions produced by unstable string-wall network would overdose the universe [38 

In the following we assume that the coherent oscillations of the saxion dominates the 
universe. This sets the initial condition for the axionic co-genesis. 

B. Dark matter and dark radiation 



The axions are produced by the saxion decay and they constitute the dark radiation. 
The axions thus produced are not thermalized because the decoupling temperature is 



about 10 9 GeV x (//10 n GeV) 



21 



40] , which is many orders of magnitude larger than 



the saxion decay temperature T a for the parameter ranges studied here. Thus the axions 



increase the effective number of neutrino species (N e g) by the amount 19|, |41] 

AiVcff = — ^ij^(^±y\ (17) 

* decouple 7 l ~ B a \g*{l*)J 

where g*(T a ) counts the relativistic degrees of freedom at the saxion decay. In deriving the 
above expression, we have used the fact that it is the entropy in the comoving volume that 
is conserved when the light degrees of freedom changes at e.g. the QCD phase transition. 
Thus, for B a ~ 0.1, 0.2, and 0.3 and g*(T a ) = 106.75, we obtain AN cS ~ 0.32, 0.71, and 
1.2, respectively. As we shall see shortly, the value of B a naturally falls in the range of 
(9(0.1), and therefore the dark radiation with AN e g = O(0.1) is a robust prediction of 
this scenario. 

Fig. [1] shows how AiV e g depends on r and the mass of the right-handed neutrino. Here 
we have assumed that the saxion decays dominantly into right-handed neutrinos, axions, 
and/or gluons, i.e., 

B a + Br(<7 -+gg) + Y, Br ^ a ^ N * N *) " L ( 18 ) 

i 

The process a ^ gg is mediated by the saxion coupling to gluons, which is analogous 
to the anomalous axion coupling to gluons (|15p . Thus Br (a — > gg) is negligibly small 
compared to B a as long as r > c^/ 71 " 2 = C(10~ 3 ). In the yellow shaded region, we 
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FIG. 1: Axion dark radiation AiV e g in the case with quasi-degenerate (left panel) and hierarchical 
(right panel) right-handed neutrinos. The amount of axion dark radiation is fixed by B a , on which 
we have imposed the condition B a + Br(cr — > gg) + J2i Br(cr — > A^iVj) ~ 1 under the assumption 
that other saxion decay modes have small branching ratios. The shaded region corresponds to 
AiV e g < 0.8, and the brown lines represent the contours for AN e Q = 0.3 and 0.6, from below, 
respectively. For the hierarchical iVj case, we have taken M\ < m CT /2 <C M2,M^. 



have AiVeff < 0.8, which is the case when B a is smaller than about 0.22: the left panel is 
obtained for quasi-degenerate right-handed neutrinos with M\ ~ Mi ~ M3 < m CT /2, while 
the right panel is for the hierarchical mass spectrum of iVj with M\ < m a /2 M 2 ,M 3 . 
The brown lines are the contours of AiV e ff = 0.3 and 0.6, respectively. 

On the other hand, the axion obtains a mass through QCD non-perturbative effects, 
and starts to oscillate coherently when the Hubble parameter becomes comparable to its 
mass. These coherent axions constitute the cold dark matter of the universe with the relic 
energy density given by 



42] 



1.19 / a n 3/4 

Q h 2 ~ 18 6 2 [ — 1 I QCD ) (19) 



where 9 a is the initial misalignment angle, and Aq CD denotes the QCD scale. Note that 
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we extracted the dependence of the QCD scale not only from the thermal effects on 
the axion mass 43| but also from the zero-temperature axion mass. Thus, the coherent 
axions produced by an initial misalignment 9 a = 0(1) make up the major part of cold 
dark matter in the universe for the axion decay constant around 5 x 10 11 GeV. For a 
larger (smaller) value of F a , we need to fine-tune the initial position close to the potential 
minimum (maximum). 



C. Non-thermal leptogenesis 

A natural way to generate a baryon asymmetry is leptogenesis via the decays of right- 
handed neutrinos. To produce a sufficient baryon asymmetry in the saxion-dominated 
universe, we consider non-thermal leptogenesis. To proceed, let us estimate the saxion 
decay temperature. This is determined by the saxion decay rate into particles that ther- 
malize to form a thermal bath: 




from which it follows, 

The above relation holds if the decay rate of Ni into leptons and Higgs is much larger 
than T a so that the right-handed neutrinos decay almost instantly after production. This 
is indeed the case when T a is smaller than m a , and the branching ratio for a — > NN is 
sizable. The decay rate of Ni into leptons is larger than r^^A^ by a factor of approxi- 
mately (/ /v) 2 m u /m a , which is about 10 2 x r (m a /T a ) 2 for B a around 0.2 and m v ~ 0.1 
eV. Thus in the parameter space with r around 0.1 and T a < m CT /2, which is of our 
interest, the relation for T a remains valid. 7 

7 More precisely, T a should be replaced by (1 — B a )T a since the axions produced in saxion decays are 
not thermalized as we will see shortly. However this changes the result only slightly for small B a . 
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The right-handed neutrinos are non-thermally produced by saxion decays when T a is 
low enough: 

777 

T, < Mi < -f, (22) 

thereby requiring T a smaller than m a /2. The lepton asymmetry generated by the decays of 
Ni is partially converted into baryon asymmetry through the sphaleron process according 
to tib = Csp^riL, where c sp h = —28/79 for T a ^> 10 2 GeV in the SM. 8 The present baryon 
asymmetry is then obtained to be 

— = c^5Zk J>Br(t7 N * N & ( 23 ) 

i 

for T a < Mi, where s is the entropy density. The condition T a < Mi is necessary to avoid 
strong wash-out of the lepton asymmetry produced by N decays. The CP asymmetry e. 



44] . and crucially 



in the N decay receives one-loop vertex and self-energy contributions 
depends on the right-handed neutrino masses and thus on /. 

Let us estimate the baryon asymmetry generated by non-thermal leptogenesis. There 
are two cases of interest, depending on the mass spectrum of Ni. First, suppose hierar- 
chical N. Then an adequate baryon asymmetry is obtained when N\ has mass not much 
smaller than m a /2 while the other two have masses larger than m a /2. This is because Aj 
is proportional to Mi, and T a less than Mi is needed to avoid strong washout effects. For 
Mi < 772,0-/2 Mi, M 3 , the CP asymmetry induced by the N\ decay is given by 



3 m Vz Mi 

ei = — 2 — ^ eff ' ( 24 ) 

l07T V 



where \S e g\ < 1 is the effective leptogenesis CP phase, and m V3 is the heaviest neutrino 
mass. Taking 5 e ff = 1, the baryon asymmetry reads 



n B 



u ( T a /m a \ / m vz x / Mi \ ( Br(a -» NiNi) \ 
10 la05evJ VlO^GevJ { 0^8 ' (25) 



The saxion mass is considered to be of order the gravitino mass. If all the SUSY particles as well as 
heavy Higgs bosons have similar masses, there are only SM particles in the plasma at the electroweak 
phase transition. 
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where an enhancement of a factor 2 can be obtained if the saxion is kinematically allowed 
to decay into the right-handed sneutrinos. On the other hand, it is also possible to have iV, 
with masses close to each other. For such quasi-degenerate iVj, the induced CP asymmetry 
is approximately given by 

(26) 



16vr v 2 AM 

where the right-handed neutrinos have masses around M with a small splitting AM. The 
CP asymmetry is enhanced by a factor M/AM, compared to the case with hierarchical 
N{. As a result, the baryon asymmetry produced in the decays of three right-handed 
neutrinos is approximately given by 



s 



10 ( T a /m a \ f m U3 \ / M \ ( M/AM \ ( B N \ 
U V 0.2 J U.05evHl0 6 GeVA 10 J \0.S ) ' 10 



for 5 cff close to 1, where B N = £\ Br(cr — > NiNi). 

The induced baryon asymmetry depends on m CT , Mi, f and B a , some of which also 
determine the amounts of axion dark matter and dark radiation. In the next subsection, 
taking into account such relations, we will explore the region of parameter space where 
the baryon asymmetry observed in the universe is explained by (|25|) or (J27J). 



D. Axionic Co-genesis 

To get a sufficient baryon asymmetry, the saxion needs to decay mainly into a pair 
of right-handed neutrinos. The relation ( fl3l) tells that this is achieved when B a is not 
small, under the natural assumption that r is not much smaller than 0.1. This implies 
that AN e g due to axion dark radiation generally lies in the range between about 0.1 and 
1 in our scenario, which is consistent with the recent Planck results (JTJ). 

Let us move on to the baryon asymmetry and dark matter. It is important to note that 
an upper bound on m a comes from the requirement that T a should be smaller than M\ for 
given /, while a lower bound arises from the requirement of generating adequate baryon 
asymmetry, provided that the right-handed neutrinos are not extremely degenerate in 
mass. The coherent axions produced by an initial misalignment 9 a = 0(1) make up 
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the major part of cold dark matter in the universe for the axion decay constant around 
5 x 10 11 GeV. For the PQ breaking scale in that range, T G will be low enough to avoid 
strong wash-out if the saxion mass is smaller than about 10 7 GeV. On the other hand, 
from (1251) and (l2Tj) . we see that ub/s ~ 10~ 10 is naturally obtained when the right-handed 
neutrino, whose decay generates the lepton asymmetry, has a mass around or larger than 
10 6 GeV. Therefore the baryon asymmetry and dark matter of the present universe can 
be achieved simultaneously for 

F a ~ 5 x 10 11 GeV and 10 6 GeV < m a < 10 7 GeV, (28) 

thereby providing a natural explanation to the baryon-dark matter coincidence problem. 
This also indicates high scale SUSY breaking around 10 6-7 GeV because the saxion is 
massless in the supersymmetric limit. Interestingly, this is the SUSY breaking scale 
needed to accommodate a 126 GeV Higgs boson within the minimal supersymmetric SM 
without large stop mixing. 

Fig. [2] illustrates how non-thermal leptogenesis works in the saxion-dominated universe. 
Here we have taken the model parameters as follows, for the cases with hierarchical and 
quasi-degenerate right-handed neutrinos: 

hierarchical N t : M x = 0Am a < M 2 , Br(cx A^iVi) = 0.9, 
quasi-degenerate jV 4 : M = 0.4m CT , M/AM = 10, B N = 0.9, (29) 

and fixed B a = 0.1 for both cases, under the assumption that r lies in the range between 
about 0.1 and 0.2. In addition, we have taken x = 0.35iVDw for the PQ charge normal- 
ization such that the PQ charges of X a are co-prime to each other, for which r is given 
by r ~ 0.12(iV DW /g) 2 . For hierarchical (quasi-degenerate) iVj with the above properties, 
ub/s = (1 ± 0.5) x 10~ 10 is obtained in the upper (lower) orange region. The observed 
baryon asymmetry can thus be explained within this region by slightly changing the prop- 
erties of Ni or 5 e ff. In the region above the gray thick line, the saxion decay temperature 
is higher than the right-handed neutrino mass, and the situation is reduced to the thermal 
leptogenesis scenario. The successful thermal leptogenesis will require the saxion to have 
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FIG. 2: Axionic co-genesis in the saxion-dominated universe. The upper (lower) orange region is 
the region where ng/s = (1±0.5) x 1CP 10 is obtained via non-thermal leptogenesis with hierarchical 
(quasi-degenerate) right-handed neutrinos. Here we have taken B a = 0.1, for which the axion dark 
radiation gives AiV e g ~ 0.32, and have assumed that the saxion decays into right-handed neutrinos 
with Bn = 0.9 (see the text). The relic axion density due to vacuum misalignment is consistent 
with the observed cold dark matter abundance in the yellow region for an initial misalignment 
angle of order unity. In the region above the gray thick line, T a is larger than M%, and the 
lepton asymmetry produced by iVj decays is strongly washed out. The dashed line represents 
T a /m a = 0.2, 0.1, 0.05, from left, respectively. 

mass of order 10 9 GeV for F a ~ 10 11 GeV. On the other hand, in the yellow region, the 
coherent axions produced around the QCD phase transition can naturally account for 
the cold dark matter in the universe. Finally we note that the axion dark radiation is 
also produced in saxion decays, for instance, AiV fj ~ 0.32 for B a = 0.1. The amount 
of baryon asymmetry does not change much for 0.1 < B a < 0.2 as long as the saxion 
dominantly decays into right-handed neutrinos. From these we conclude that the axonic 
co-genesis is realized in the region where the orange and yellow bands meet, that is, for 
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F a around 5 x 10 11 GeV and m a in the range 10 6-7 GeV. 

Before closing this section, let us comment on the lightest supersymmetric particle 
(LSP) which is stable under the R-parity conservation. In high scale SUSY, the thermal 
relic abundance of the LSPs may be too large. This can be avoided if T a is lower than 
the LSP freeze-out temperature. If LSPs are non-thermally produced by saxion decays, 
its branching fraction must be suppressed. Another way is to consider some mediation 
mechanism rendering the LSP sufficiently light, or to introduce a very light superparticle 
in a hidden sector into which the ordinary superparticles decay. Alternatively one can 
simply assume R-parity violation. We also note that the results on the baryon asymmetry 
and dark matter in the saxion-dominated universe hold even in non-supersymmetric axion 
models, if the saxion is stabilized with m a <C /. 



IV. SAXION STABILIZATION 



In this section we discuss how to make the saxion dominantly decay into right-handed 
neutrinos rather than into axions. This requires a small r in ( 1131) . If PQ symmetry is 
spontaneously broken by a single field X, which is the case for instance in the radiative 
saxion stabilization, one obtains r = 1, implying that Br(a — > iVjiVj) is much smaller than 
B a . Even in this case, one can avoid an axion-dominated universe by introducing the 
superpotential term, X 2 H u Hd, so that the saxion dominantly decays into SM particles 
or higgsinos. However small Br(cr — > NiNi) makes it difficult to implement non-thermal 
leptogenesis successfully. This is resolved when multiple fields are involved in the PQ 
symmetry breakdown. 

Let us consider supergravity models where two PQ fields, X\ and X 2 , obtain VEVs 
around /. One simple way to stabilize the saxion is to consider 

W = y x Y.(X x X 2 -n% (30) 

with qi + q 2 = and \i ~ /, where £ is a PQ singlet field. In this case, the F-flat direction 
X\X 2 = fi 2 is lifted by soft SUSY breaking terms, 

-C aait = ml\X 1 \ 2 + m*\X 2 \ 2 . (31) 
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For yx^ 3> rn ?J / 2 and positive m 2 2 similar to or smaller than rn. 2 ^, the PQ fields are 
stabilized at 

m\v{ ~ m\v\, (32) 

with V1V2 = yU 2 ! . Thus, if m 2 = m|, the PQ fields are fixed with v\ ~ for which r 
has a tiny value. This implies that r can have a small value when X% and X 2 feel SUSY 
breaking via couplings of similar strength. For Xl j2 having soft scalar masses m 2 ± Am 2 
with Am 2 <C m 2 , we get r ~ Am 2 /m 2 . If Am 2 /m 2 = 0.1, which would not require severe 
fine-tuning of model parameters, the relation ( fl3|l tells Br(a — > NiNi) < 3.6B a . In this 
case, one can arrange for instance Br(er — > NiNi) around 0.7 and B a ~ 0.2 by taking 
appropriate right-handed neutrino masses. Note that + B a is not necessarily equal to 
one in the presence of a superpotential term, X 2 H u H d . 

Another way to stabilize the saxion is through the competition between SUSY breaking 
effects and a higher dimensional superpotential term: 

W = k f\ X \ , (33) 

with n\ + n 2 > 4 and q\U\ + g 2 n 2 = 0. Here we take q% + g 2 7^ to forbid a large mass 
term in the superpotential. It is straightforward to find 

= -q 2 v 2 Jl + OC^X) ~ (A- + -- 3 m 3/2 ) 2/(ni+n2 - 2) , (34) 

V Vm 3/2 

and m a ~ m 3 / 2 , assuming |m 2 2 | <C ?7ig/ 2 . The value of r is fixed to be r ~ (1 — ni/n2) 2 
or (1 — n 2 /n!) 2 , depending on which PQ field is responsible for Mj. Let us suppose that 
iVj obtain masses from the coupling to one of the PQ fields that gives a smaller value of 
r. Then the saxion stabilization by the higher dimensional operator leads to 

(m,n 2 ) = (1,3) : r~0.44, / ~ A 1/2 m^ 2 , 

(ni,n 2 ) = (2,3) : r ~ 0.11, / ~ A 2/3 m^. (35) 

The first model has Br(a — > NN) < 0.8B a and B^ < 2AB a . Thus, to get a sufficient 
baryon asymmetry, one would need heavy N\ with mass close to 10 7 GeV if the right- 
handed neutrino masses are hierarchical. One may instead consider degenerate Ni with 
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M ~ 10 6 GeV and M/AM ~ 10, for which the successful non-thermal leptogenesis is 
achievable with B a ~ 0.3 and Bn ~ 0.7. On the other hand, in the second model, we have 
Br(<j — > NiNi) < 3.6B a and Bn < 10.8i? a . It is thus easy to make the saxion dominantly 
decay into right-handed neutrinos. We also note that, in both cases, the PQ breaking 
scale / naturally lies in the range 10 n ~ 13 GeV for 777,3/2 ~ 10 6-7 GeV, taking A to be the 
GUT or Planck scale. 9 

Our discussion so far has assumed that the saxion given by (jl]) does not mix with other 
CP-even scalar bosons, and that its decay rate into the PQ sector scalars or fermions is 
small compared to that into right-handed neutrinos. For the model where X\ and X2 
are stabilized by a higher dimensional superpotential term, there appear additional light 
scalars and singlino. Hence more care should be taken. The saxion a indeed corresponds 
to the eigenstate in the limit of vanishing soft scalar masses, and there is induced small 
mixing with the other CP-even scalar for \m\ 2 \ ^3/2- Assuming small soft scalar 
masses, we find that the additional CP even and odd scalars are both heavier than the 
saxion in the models ( 13 5p . The PQ sector also includes two fermions composed of X\ and 
Xi- They have masses about 2.Ylm a and 0.71m<j for (711,712) = (1)3), while 2.83m CT and 
0.71m CT for the model with (711,712) = (2,3). Thus the saxion decay into PQ fermions is 
kinematically closed in both models. 



V. DISCUSSION AND CONCLUSIONS 



The PQ symmetry may be unbroken during inflation. Then the spontaneous break- 
down of the PQ symmetry after inflation leads to the production of topological defects 
such as axionic strings and domain walls. According to the recent numerical simulations 



of axion string- wall network 38|, |39[, the PQ breaking scale should be of order 10 10 GeV 



for explaining the dark matter. Based on our above argument, it will be difficult to 



9 The model discussed in Ref. Z7\ also realizes an intermediate axion decay constant by relating it to 
the SUSY breaking scale, F a ~ yJm^J^Mpl. However it is difficult to implement axionic co-genesis in 
this model because the PQ sector itself participates in SUSY breaking, leading to a too large saxion 
mass around F a . 



18 



combine this scenario with non-thermal leptogenesis by the saxion decay. See Fig. |2j 
However, thermal leptogenesis will be possible, if the reheating temperature is sufficiently 
high Tr > F a ~ 0(1O 10 ) GeV. 10 The axion will be thermalized, but its contribution to the 
effective number of neutrinos is suppressed because of large relativistic degrees of freedom 
at high temperature. Note that the saxion does not have to dominate the universe in 
this case, and the application to a non-SUSY case is straightforward. It is interesting 
to note that the baryon and dark matter abundances are determined by a single scale 
Mi ~ F a ~ O(10 10 ) GeV. 

If the PQ symmetry is broken during inflation, the axion acquires quantum fluctuations, 
generating isocurvature perturbations. If the PQ breaking scale F a does not change 
during and after inflation, the inflation scale is constrained by the isocurvature constraint 
as i?i n f < (9(10 7 ) GeV. This bound is relaxed if one of the PQ scalars takes a larger 
field value during the inflation. Such large deviation will ease the conditions for the 
saxion domination. The isocurvature perturbations of the axion dark matter may be 
found or constrained by the future observations of the CMB power spectrum and non- 



Gaussianity [46l-l51|. 



In this paper we have considered the PQ sector as the origin of baryon asymmetry, 
dark matter and dark radiation; the axionic co-genesis. We have found that non-thermal 
leptogenesis by the saxion decay works successfully for the saxion mass of O(10 6 ~ 7 ) GeV 
and the PQ breaking scale around 5 x 10 11 GeV. The dark matter can be explained by 
the axions produced in the misalignment mechanism with an initial deviation of order 
unity. If the existence of dark radiation is confirmed by future observations such as the 
Planck satellite, this scenario will be one of the plausible solutions to the cosmological 
coincidence problems. Intriguingly, in the SUSY framework, the suggested SUSY breaking 
scale is consistent with the observed SM-like Higgs boson of mass 126 GeV in the minimal 
supersymmetric SM without large stop mixing. Thus, the axionic co-genesis may give an 
answer to the question of why the SUSY breaking scale is much higher than the weak 



10 



The thermal leptogenesis in the axion (=majoron) model was considered in Ref. [fj. See also Ref. [4£ 
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scale, which generically requires severe fine-tuning for the correct electroweak symmetry 
breaking. 
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